function [ s ] = evaluateQuadraturePoint( p, q, qn )
%EVALUATEQUADRATUREPOINT Summary of this function goes here
%   Detailed explanation goes here

r = p - q;

rn = norm(r);

rg = r / rn;

% qnc = [qn(2); -qn(1)];

M = [rg(1)*rg(1), rg(1)*rg(2), rg(1)*rg(3); rg(2)*rg(1), rg(2)*rg(2), rg(2)*rg(3); rg(3)*rg(1), rg(3)*rg(2), rg(3)*rg(3)];

chi = cross(qn, rg);
N = [0, -chi(3), chi(2); chi(3), 0, -chi(1); -chi(2), chi(1), 0];

I = diag(ones(3, 1));

prpn = dot(rg, qn);

s = -1 / (8 * pi * rn^2) * (prpn * (I + 3 * M) - N);

end

